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Abstract — In this paper, we investigate the minimum average 
transmit energy that can be obtained in multiple antenna broad- 
cast systems with channel inversion technique. The achievable 
gain can be significantly higher than the conventional gains 
that are mentioned in methods tike [1]. In order to obtain this 
gain, we introduce a Selective Mapping (SLM) technique (based 
on random coding arguments). We propose to implement the 
SLM method by using nested lattice codes in a trellis precoding 
framework. 

I. Introduction 

Recently, there has been a considerable interest in Multi- 
Input Multi-Output (MIMO) antenna systems due to achieving 
a very high capacity as compared to single-antenna systems. 
Multiuser MIMO systems can also exploit most of the advan- 
tages of multiple-antenna systems. 

In a broadcast system, when an access point with multiple 
antennas is used to communicate with many users, the commu- 
nication is complicated by the fact that each user must decode 
its signal independently from the others. As a simple precoding 
scheme, the channel inversion technique can be used at the 
transmitter to separate the data for different users. However, 
this method is vulnerable to the poor channel conditions. 

In this paper, we investigate the optimum gain for average 
transmit energy in multiple antenna broadcast systems with 
channel inversion technique. By using the fact that the channel 
is not orthogonal, the gain that can be achieved is significantly 
higher than the regular shaping gains that can be achieved in 
methods like [1]. 

In a broadcast system with the channel inversion tech- 
nique (given a fixed channel matrix), we find the optimal 
probability distribution for the data vectors to minimize the 
average transmit energy. Then, we introduce a theoretical 
Selective Mapping (SLM) technique (based on random coding 
arguments) to obtain the optimal average transmit energy. In 
order to implement the SLM method effectively, we propose 
using nested lattice codes in a trellis precoding framework 

The rest of the paper is organized as follows. In Section II, 
the system model is introduced. Section III finds the optimal 
probability distribution for transmit data in channel inversion 
techniques. Section IV is devoted to introducing the SLM 
technique and its analysis and implementation issues. 
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II. System Model 
A multiple antenna broadcast system can be modeled by [2] 

y = Hx + n, (1) 

where y is the Mx 1 received vector, x is the Nx 1 normalized 
transmitted data,n is additive white Gaussian noise, and H 
represents the M x N channel matrix in real space. 

In broadcast systems, the receivers should decode their 
respective data independently and without any cooperation 
with each other. The simplest method is using the channel 
inversion technique as a precoding method at the transmitter 
to separate the data for different users s = H + u, where 
H+ = H*(HH*) \ H* is the Hermitian of H, u is the 
data vector, i.e. Ui is the data for the i'th user, and s is 
the transmitted signal before the normalization. When the 
number of transmit antennas equals with the number of users, 
M = N := M, the transmitted signal is 

s = H ^u. (2) 

As in [1], the normalized transmitted signal would be x = 

g 

— , where 7 = ||s|| 2 . The problem arises when H is 

poorly conditioned and 7 becomes very large, resulting in a 
high power consumption. 

In a multiple antenna system, it is assumed that the data 
vector u is selected from a constellation with discrete points. 
However, through this paper we investigate the probabilis- 
tic behavior of the transmitted signal s. Assuming a large 
constellation, continuous approximation provides a probability 
distribution for each constellation, resulting in different E{^}. 
The challenge is finding the best probability distribution with 
minimum E{j}. Note that the expectation in E{^f} is over u 
and the channel is assumed constant. 

III. Optimum Probabilistic Constellation 

Channel inversion technique removes the need for complex 
decoding algorithms in the receiver side; however, it leads to a 
high energy consumption as the average energy of the resulting 
constellation points is high. We are looking for a constellation 
shaping method for the input constellation, such that using the 
channel inversion technique, the resulting constellation has a 
smaller value for the average transmit energy. 



A proper input constellation should be designed such that 
two conditions are satisfied: (i) data can be decoded indepen- 
dently at the receivers, and (ii) the average transmitting energy 
is as low as possible. The design of the constellation is known 
as the shaping. By using a conventional block constellation, 
any point in the constellation is equally likely. However, 
by shaping, a nonuniform distribution is achieved over each 
dimension. 

A common constellation shaping technique is to choose a 
finite set of points from an M dimensional lattice A that lies 
within a finite region TZ C R M . This constellation is known 
as a lattice code. If C is a lattice code of reasonably large size, 
then the distribution of its points in M dimensional space is 
well approximated by a uniform continuous distribution over 
the region TZ (the continuous approximation) [3]. 

Having a uniform distribution over region TZ induces a 
nonuniform distribution on each dimension. In other words, 
if u is selected uniformly over 1Z, each element of u has 
a nonuniform distribution. Through this paper, the probability 
distribution of the elements of u is called marginal probability 
distribution of u. We assume that the region 1Z has a fixed 
volume Vol(7?.) = V, resulting in the entropy of (logV). In 
the case of independent variables, we assume that the entropy 
per real dimension is Ti — log V. 

Let Q := (H _1 ) T H _1 = UAU T , where U is the unitary 
matrix of eigenvectors of Q and A is the diagonal matrix of 
the corresponding eigenvalues, Xi, i — 1, ••• , M, Assume 
u e 1Z be a random vector with mean ^{u} = /x and the 
correlation matrix E{uu T } = X > 0. The energy of the 
transmitted signal, called transmit energy, is defined as 7 = 
u T Qu = u T UAU T u = v T v, where v = \/AU T u. When u 
is selected uniformly in TZ, the vector v is selected uniformly 
over a region 1Z', where TZ' = {v|v = \/AU T u,Vu 6 TZ}. It 
is more convenient to explain some behaviors of 7 based on 
v. The average transmit energy can be written as [2] 

£{7} = tr(QE) + M T Q/x (3) 

If we ignore that users are supposed to decode their data 
independent of each other, the optimum region for the input 
constellation can be found using the following lemma: 

Theorem 1: Let u — [ui,U2, ■ ■ ■ , mm] S K- M be a random 
vector with probability distribution f(ui, 112, ■ ■ ■ , %/), mean 
E{u} = fi, and the correlation matrix £{uu T } = X > 0, 
in a broadcast system introduced in (fTJ. Let H(u) denote the 
entropy of the data vector u. Then, a multivariate Gaussian 
random vector u with fi — and the covariance matrix 

s = \^nV7 2 HH T (4) 

will minimize the energy of the transmit signal given a fixed 
entropy H(u) = log(V), where a 2 is the variance of a 
Gaussian random variable with entropy H = jj log(V). 

Proof: See [2] ■ 
This choice of S suggests that the minimum value of the av- 
erage energy among transmit signals with different probability 
distributions is [2] 

£{7} = M h {/u\a 2 (5) 



Consider the auxiliary vector v = VAUu. It can be 
easily shown that each element of this vector has a Gaussian 
distribution with zero mean with variance TZ 2 ,. Therefore, 
in the limit of M — ► 00, this vector is uniformly selected 
over an M-dimensional sphere centered at the origin with 



radius v AITZ eq , i.e. B m (0 , V MlZ eq ) (corresponding to the 
minimum average transmit energy in (0. 

Roughly speaking, we can assume that the region TZ' is 
(0, ^/MTZ eq ). On the other hand, the vector u, a Gaussian 
random vector with zero mean and covariance matrix in (|4), 
is uniformly selected over the region TZ which is an oval. The 
main diameters of this oval are along the eigenvectors U and 
the radii of the oval in each direction are \J~^-TZ eq for i = 

1, • • • , M, in other words TZ = O M (0, ^j^TZ eq ) . (see [2]). 
By using this region, an additional channel gain of [2] 

Arithmetic Mean(Xi, • • • , Ajv) 



Gh 



(6) 



E, 



opt 



Geometric Mean[X\, ■ ■ ■ , Ajv) 

can be achieved (in addition to the conventional shaping gain). 

The geometric mean of a data set is always smaller than or 
equal to the set's arithmetic mean (the two means are equal if 
and only if all members of the data set are equal). On the other 
hand, without the channel matrix, we have the conventional 
shaping gain. However, the presence of H 1 will affect the 
shaping gain by the Channel Gain, Qu, defined in (|6). Without 
the channel effect the optimum region TZ is a spherical region 
(corresponding to independent Gaussian variables), while with 
the channel effect the optimum region TZ is an A/-dimensional 
oval. 

From another point of view, this gain can be seen as the 
effect of rate (or power) allocation for Gaussian distribution 
which has been considered in multi-carrier transmission and 
point to point multiple antenna systems, e.g. [4]. However, 
this concept ignores the independency condition required for 
a broadcast system. Here, the challenging problem is how 
the region TZ or TZ' can be achieved, while considering the 
independency condition. 

IV. Selective Mapping 

The idea of Selective Mapping (SLM) is to generate a large 
set of data vectors that represent the same information, where 
the data vector resulting in the lowest energy is selected for 
transmission. This idea has been used in OFDM systems, e.g. 
[5], to reduce the average transmit energy. 

In the optimum case, the vector u is selected uniformly 
over an M-dimensional oval and the transmit vector is selected 
uniformly over a hypersphere. However, due to the indepen- 
dency condition, implementing this oval shape region is not 
possible. The receivers can not co-operate with each other to 
locate a point inside this oval. We propose an SLM method 
that can theoretically achieve the optimum gain for average 
transmit energy. The region for vector u is not oval; however, 
the resulting region for the transmit vector in the limit is a 
hypersphere. 

In the sequel, first, we use a random coding argument 
to explain the SLM method, its analysis, and the maximum 



theoretical gain that can be achieved. In this part, again we 
ignore the independency condition. In continue, we implement 
the SLM technique considering the independency condition by 
using a trellis precoding. 

In the system model (Q]i, the volume of the region is fixed, 
V. In order to provide multiple choices for the SLM method, 

the volume is increased to V such that for each data vector 

V 

there are N points, where N = — . 

V 

In other words, N i.i.d. samples of u are generated, 
{ui, U2, • • • , Ujv}, and s; with the lowest transmit en- 
ergy is selected for transmission. In other words, ji = 
min{7i, 72, • • • , 7jv}- We are looking for the probabilistic 
behavior of 72 . 

A. Asymptotic Analysis 

In this section, we analysis the effect of an SLM method 
for broadcast systems. In the proposed method, N i.i.d. 
samples of u are generated, {ui,U2,--- , u/v}, and among 
the corresponding transmit vectors Sj = H x Ui, the vector s; 
with the lowest transmit energy is selected for transmission. 
In other words, in the SLM method, we are looking for 



mm Sj 

Ki<N 



(7) 



where ||.|| represents the regular norm. 

The expression in (0 is similar to minimization of distor- 
tion in quantization and random quantizers. The tremendous 
research on random quantization [6, and ref. therein] can help 
us to evaluate the expression in (|7]i in our SLM method. 

Let si,S2,-- - ,sjv, be i.i.d. R M -valued random variables 
with distribution Q, i.e. 

Q(v)=P{ Sil <«!,.•• ,s iM <v M } i = l,---,N, (8) 



where 



(vi, ■ ■ ■ ,v M ) e m 



A I 



For any region 1Z, the probability Q(1Z) is the probability that 
there is at least one code point in the region 1Z, i.e. 



Q(K) 



Q(dy). 



n 



Define the r th order transmit energy as 



7 r ,iV 



mm s, 

Ki<N 



7a> 
Q 



(9) 



In this 



where based on our previous notation 7; 
section, the asymptotic probabilistic behavior of J2N> wnen 
N — ► 00, is investigated. Specifically, we calculate the 
average transmit energy in the SLM technique. Note that, 
in the following, we frequently use A which is the M- 
dimensional Lebesgue measure. Here, we define it as the M- 
dimensional volume of a region. 

Theorem 2: Let 81,82,- •• , Sjv, be i.i.d. R M -valued ran- 
dom variables with distribution Q. Then, 



lim E 



= B 



M 



T(l + —)a p 



(10) 



where B x = 2, B M = A (B M (0, 1)) = tt m / 2 /T(1 + M/2) for 
M = 2, • • • , and g p is defined for any p > as 

Q(B M (0,S)) 



9p 



inf 



ae(o, P ] X(B M (0,6)) 
Proof: See [2]. ■ 
Now, consider the special case of uniform distribution. 
When we have a large lattice code, we can assume we have a 
uniform distribution over the region where the lattice code is 
defined. Applying SLM technique, over a region with uniform 
distribution results in the following average for the r th order 
transmit energy. 

Theorem 3: Let 1Z C R M be a compact set with X(TZ) > 
and let si,-- - ,sn be i.i.d. random variables with uniform 
distribution over 1Z. Then, 



lim E 

JV— 00 



= B M M T(l + —)X(n)^. (11) 

Proof: Let Q be a uniform distribution over 1Z, i.e. Q = 
U{1Z). Therefore, 



B M (0,i-r))=- 



X(K) 



and 



9 P 



= inf 



<(B M (0,S)) 



1 



(12) 



(13) 



sg(o, p ] A(B M (0,5)) X(1Z) 

Substituting ( fT3l in ( TTOb completes the proof. ■ 
Note that we are interested in cases that the i.i.d random 
variables ui, ■ • • , un are selected uniformly over a region 1Z'. 
According to s = H _1 u, for the probability distribution of s, 
we have 

/ s (s) = IH- 1 ! / u (Hs). (14) 

Therefore, if u has a uniform distribution over 1Z', s has also 
a uniform distribution over K, where 1Z = H^ 1 ^'. 

In order to find the asymptotic average transmit energy of 
SLM technique, we should replace r — 2 and Q — U(1Z), 
where 1Z is the region for the transmit vector s. Therefore, 
according to the expression in ( fTTT i. the average transmit 
energy for large N can be approximated by [2] 

Eslm =T(1 + ^)MTl 2 eq . (15) 



M 



Comparing <(5j and (fTBI l, we can see that using SLM technique 
with any lattice code of reasonably large size the optimum 
transmit energy can be achieved since for large M, T(l + 

&) — r(i) = 1. 

Corollary 1: In a broadcast system, applying SLM method 
to lattice codes of reasonably large size, with a fixed volume, 
will result in equal values for the average transmit energy when 
N is large enough in the SLM method. 

We must emphasis that in our random coding argument the 
probability of the event that two different code words have 
the same transmit data vector is negligible. In the case of this 
event, we have an error in our broadcast system. However, 
since the probability of this event is small, the average transmit 
energy would not change. 



B. Implementation Issues 

In any practical SLM method, the lattice code C (constella- 
tion) should be expanded such that the number of constellation 
points are multiplied by N, resulting in a new lattice code 

C. This new set of constellation points are grouped in |C| 
sets containing N points. Transmitting any of these N points 
transfer the same information. These sets (and the expanded 
constellation) should be selected such that the users at the 
receive side can decode their data independent of each other. 

The method proposed in [1] can be considered as an SLM 
technique with this idea. In this method the region for the 
transmit vector s is expanded by repetition of the constellation 
by multiples of r in each direction. In other words, for any 
vector s = H V we find H _1 (u+rl) for [-b/2\ +1 < k < 
[b/2\. In each direction we repeat the constellation b times, 
so we have N = b M in the SLM method. In this method, the 
transmit vector s is selected in the original constellation and 
N — 1 other points are calculated by adding integer vector 
offsets, resulting in N points in the expanded lattice code C. 
A modulo operation in the transmitter and receivers guarantees 
and satisfies the independency condition. 

For large enough M and N, this method can not achieve the 
optimum average energy. The equivalent region for vector s is 
the Voronoi region of tH -1 , not a hypersphere [2]. This leads 
to an improvement over the Gaussian marginal probability 
distribution; however, this is not the best that we can achieve. 
The more this Voronoi region looks like a sphere, the less the 
average transmit energy is. The problem in the SLM method in 
[1] is that the vector s is not uniformly distributed over lattice 
code C. In order to preserve the independency condition, a 
vector is uniformly distributed over C and N — 1 other points 
are calculated deterministically in C based on this point. This 
results in a region with Voronoi region shape not a sphere. 

In [7], a sign-bit shaping algorithm is proposed for pre- 
coding in broadcast systems. Sign-bit shaping is implemented 
by using a trellis code. This technique is actually an SLM 
method since it gives the transmitter many different options 
when determining which symbol to transmit. 

Trellis shaping systems are composed of a rate (k s ,n s ) 
binary convolutional shaping code C and a signal set A 
partitioned into 2 n " shaping subsets [8]. The signal set A is 
typically a lattice code with shaping region 1Z, and the shaping 
subsets are the points of this region that fall within subregions 
IZi for i = 1, • • ■ , 2" a . It is important that the code C and the 
signal set A is selected such that the the equivalent points are 
selected uniformly over 1Z. 

Conventional precoding schemes in broadcast systems, such 
as [1], treat multiple antennas of different users as different 
users. By using trellis shaping for each virtual user, in each 
2-dimensional space, there is a modulo operation with respect 
to the Voronoi region of the shaping trellis code. In other 
words, in the space of each user, there is a modulo operation 
with respect to the Cartesian product of these Voronoi Regions. 
However, we can use the shaping concept in each user's space. 



The idea in [7] can be extended to include the multi-antenna 
case. We can use a trellis shaping for each user, and not for 
each antenna. 

Nested lattice codes can be implemented such that both 
these improvements are met. The idea of nested lattice codes 
has already been used for interference cancelation in de- 
generated broadcast systems [9]. There, it is assumed that, 
in an ordered set of users, each user has the ability that 
it can decode the message for the previous users. In other 
words, it is assumed that each user has the code-book for 
the previous users. This technique can be implemented in our 
scheme to provide us with shaping, without any need for these 
assumptions. We can achieve the same gain as that reported 
in [9] for broadcast systems with precoding, without any extra 
assumption. 

Assume that in a broadcast system with K users, each 
user has n u antenna in (Q]i, i.e. M — 2Kn u . One way 
of implementing this idea is implementing a large trellis 
consisting of K sub-trellises with a lattice partition, A/ A', 
in a 2n u dimensional space. In each sub-trellis, the lattice 
code is divided into |A/A'| partitions, i.e. for transmitting any 
information vector, one of IA/A'1^ equivalent points, in M 
dimensional space, with the lowest transmit energy is selected. 
In other words, the vector u resulting in the lowest energy for 
H _1 u would be selected for transmission. Now, in each 2n u 
dimensional space, the modulo operation is with respect to the 
Voronoi region of this trellis code. The Cartesian product of 
these regions should be as close as possible to the sphere in 
order to generate the optimum shaping region. 
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